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ABSTRACT 


This  paper 


C..  • 


f-'a  "serte^B  of  three  which  analyze  an  adaptive  approxi¬ 


mate  approach  for  solving  (n+l)-dimensional  boundary  value  problems  by  replacing 

them  with  systems  of  equations  In  n-dlmenslonal  space. 

'TXt  <r7^y 

In  this  'peper-wd’  show  how  to  find  reliable  a-posterlorl  estimates  for  the 
error  and  how  these  can  also  be  used  in  the  design  of  an  adaptive  strategy. 


Various  numerical  examples  are 
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INTROTUCTION 


In  a  recent  paper,  [6],  we  Introduced  the  concept  of  dimensionally  reduced 
solutions  to  an  elliptic  boundary  value  problem.  These  are  obtained  by  project¬ 
ing  (In  the  energy)  the  true  solution  of  the  boundary  value  problem  in  the  n-fl- 
dlmenslonal  domain  u  x  [-h,hl  onto  spaces  of  the  form 

h  ** 

^  1  w-Cx))!)  (y/h)|w  arbitrary}  , 

j»0  J  J  ^ 

00 

where  'Is  a  given  set  of  functions  on  [-1,1],  (x  are  coordinates  on 

and  y  ranges  over  [-h,h]).  For  some  basic  Ideas  behind  this  concept,  see 
[6]  and  the  Introduction  to  [5].  In  [6]  the  focus  was  on  the  right  selection 
of  the  ij^j's  .  It  was  shown  there  that  for  a  very  wide  class  of  problems  the 
'|»j's  should  be  selected  such  that 

span  , 


where  P  Is  a  second  order  ordinary  differential  operator  intrinsic  to  the 
elliptic  boundary  value  problem. 

In  [7]  we  analyzed  the  convergence  properties  of  such  methods  as  the  order, 

N  ,  increases. 

The  present  paper,  which  Is  a  direct  continuation  of  the  previous  work, 
deals  with  the  problem  of  reliable  a-posteriori  error  estimation.  It  also  de¬ 
signs  an  adaptive  algorithm  for  the  selection  of  the  right  dimensionally  reduced 
solution.  As  It  follows  from  [6]  and  [7]  a  high  number  of  basis  functions 
may  be  needed  (depending  on  the  desired  accuracy)  either  if  the  thickness  of  the 
domain,  h  ,  Is  not  sufficiently  small  or  there  are  singularities  In  the  true 
solution  to  the  boundary  value  problem.  (Because  of  the  corner  In  the  domain  such 


2 

singularities  are  often  present  in  the  neighborhood  of  Sw  {-h}  and 
B'o)  {h} . ) 

Since  singularities  are  local  phenomena  it  is  of  utmost  practical  Importance 
to  introduce  dimensionally  reduced  solutions  that  permit  N  ,  the  order,  to 
vary  throughout  the  domain  u  .  This  aspect,  specifically  the  adaptive  choice 
of  the  distribution  for  N  ,  is  also  addressed  here. 

We  now  give  a  short  review  of  the  contents  of  this  paper. 

In  Section  2  we  give  a  precise  formulation  of  the  model  problem  (which  is 
identical  to  that  of  [6])  and  prove  some  auxiliary  results. 

Section  3  is  devoted  to  the  construction  of  an  estimator  for  the  error. 

The  main  theoretical  results  in  this  section  are  Theorem  3.1  and  Theorem  3.2, 
which  show  that  the  introduced  estimator  is  an  upper  bound  for  the  error  but  on 
the  other  hand  is  not  too  conservative  (  away  from  singularities  and  for  rea¬ 
sonably  small  h  ) .  Numerical  experiments  verify  this  and  furthermore  indicate 
that  even  for  relatively  large  h  or  strong  singularities  the  estimator  is  of 
the  same  magnitude  as  the  error.  The  problem  of  how  to  detect  if  the  estimator 
is  unacceptably  conservative  and  how  to  improve  it  are  addressed  in  Section  4. 

In  Section  5  we  extend  the  concept  of  dimensional  reduction  to  include  a 
possibly  different  number  of  basis  functions,  ,  in  different  parts  of  the 
domain.  We  also  design  an  adaptive  strategy  to  select  the  right  distribution  for 
the  number  of  basis  functions.  This  strategy  is  based  on  our  ability  to  give 
reliable  estimates  for  the  error  much  in  the  same  way  as  the  strategy  used  by  the 
finite  element  solver  F.E.A.R.S.  to  generate  an  'optimal'  grid  (cf.  [1]). 

Finally  Section  6  (and  also  4)  contains  a  numerical  example  that  Illustrates 
how  well  the  error  estimation  and  the  adaptive  strategy  perform  in  practice. 
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2.  NOTATION  AND  THE  MODEL  PROBLEM 

Let  H  be  a  separable  Hilbert  space  with  Inner  product  <u,v>  and  norm 
Mull  “  <u,u>  . 

A  denotes  a  (possibly  unbounded)  self-adjoint  linear  operator  In  H  with 
domain  of  definition  D(A)  . 

Furthermore  we  assume  that  A  Is  a  strictly  positive-definite  operator, 
i.e.,  there  exists  C  >  0  such  that 

Vu€D(A):  C|  [uj  _<  <Au,u> 

Let  M  be  a  self-adjoint  bounded  linear  operator  in  H  .  M  is  also 

assumed  to  be  a  strictly  positive-definite  operator. 

1/2  1/2  1/2 
V(A  )  Is  Itself  a  Hilbert  space  with  inner  product  <u,v>  +  <A  u,A  v> 

-1  k 

The  same  is  true  about  P((M  A)  )  for  any  Integer  k  ^  0  . 

2 

I  denotes  an  interval  on  the  real  line.  L  (I;H)  is  defined  as  the  set 

of  strongly  measurable  functions  u:  I  H  such  that  ||u(*)||  is  an  element 

of  L^(I),  (cf  [4]).  The  same  goes  for  L^{I;I>(A^^^))  and  L^(I;P((M“^A)*^))  . 

1/2 

We  also  need  Sobolev  spaces  of  functions  with  values  in  H,P(A  )  and 
—  1  k  1 

P((M  A)  )  .  H  (I;H)  denotes  the  space  of  functions  u:  I  -►  H  such  that 

u(')€L^(I;H)  and  u(’)€L^(I;W)  ,  (cf.  [2]).  The  spaces  for  P(A^^^) 

—1  k 

and  P((M  A)  )  are  defined  similarly.  The  derivative  is  taken  in  the  distri¬ 
butional  sense. 

H^(I)  denotes  the  standard  Sobolev  space  on  I  . 

go 

Assume  a  and  b  are  real  valued  functions  in  L  ([-1,1])  such  that 


Uq  - 

bg  1  b(y) 


4 


00 

for  some  constants  a  >  0  ,  b  >  0  .  a.  and  b.  €L  (t-h,hl)  are  then  defined 

o  o  n  n 


aj^(y)  -  a(y/h) 
bj^(y)  *  b(y/h) 


By  denote  the  differential  operator  -  ■^)  •  Let  f  and 

g  be  two  arbitrary  vectors  from  H  .  We  consider  the  following  model  problem 


0 

in 

l-h,hl 

d  „  h 

\d7”“  • 

g 

for 

y  *  h 

d  „  h 

%d7”“  ■* 

f 

for 

y  »  -1 

The  precise  formulation  of  (1)  is 


u^€H^([-h,h];H)nL2([-h,h];!?(A^''^)) 


8^(u  ,v)  =  <g,v(h)>  -  <f,v(-h)> 


.V'v€H^([-h,h];H)nL2([-h,h];P(A^/^)) 


where  8,  denotes  the  bilinear  form 
n 
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For  more  details  see,  [6}.  In  that  paper  we  Introduced  the  notion  of  dimensionally 
reduced  solutions  to  (1).  Let  ([-1,1])  be  a  given  sequence  of 

linearly  Independent  functions  (referred  to  as  basis  functions) . 

Definition 

The  dimensionally  reduced  solution  ujj  of  order  N  Is  the  projection 
of  u^  onto  the  space 

h  ^  wo 

V  =  {  J  i()  (y/h)xjx,  €P(A^'‘^),  j  =  0,..,N} 

j=0  J  J  J 

The  projection  is  with  respect  to  the  inner  product  Bj^(u,v)  . 


We  proved  that  in  order  to  obtain  optimal  rate  error  estimates  for  h  0 


there  Is  essentially  only  one  choice  for  the  sequence  •  This  is  related 

to  the  operator  P  =  b  ^  ^  a  ^ 


THEOREM.  There  exists  a  sequence  of  linearly  Independent  functions 
.  with 

(I)  W(P^)  =  span{i()j  1  >  1  , 

that  has  the  following  property: 

(II)  For  any  Integer  N  ^  0  and  for  any  given  set  of  vectors 

f.g€P((AM  ^)”)  there  exists  a  constant  (Independent  of  h) 
such  that 


h  h  , , I  „  ,  2N+1/2 


11^  -“2nII!e  V 


W(P^)  here  denotes  the  nullspace  of  P^  ,  and  |||•||L  Is  the  energy-norm 
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associated  with  the  bilinear  form  8j^  .  This  is  slightly  different  from  the 
formulation  in  [6],  where  we  used  the  norm 

(|  1  1-^  u(y)  I  +  I  I  lA^'^^u(y)n^dy)^'^^ 

-h  -h 

It  Is  obvious  though,  that  these  two  norms  are  equivalent  with  constants  in¬ 
dependent  of  h  . 

For  more  details  concerning  this  theorem  and  its  converse  we  refer  to  [6]. 
It  is  now  conveinent  to  Introduce 


Definition 

-Oo 

Any  sequence  *^^*‘*-  properties  listed  in  the  previous 

theorem  is  said  to  be  an  optimal  sequence  of  basis  functions. 


It  follows  immediately  from  Theorem  4.1  of  16]  that  any  two  optimal  sequences 

oo  oo 

of  basis  functions  {♦,}.  «  and  {i|/. },  ~  satisfy 

J  j“0  J-0  ^ 


span{4ij}j^Q  -  span{iJij}”^Q  VN  1  0  . 


We  shall  often  use  this  fact  without  explicitly  mentioning  so. 

In  the  present  paper  we  need  a  slightly  different  but  weaker  version  of  the 
result  contained  in  Theorem  4.1  of  [6]. 


LEMMA  2.1.  Let  optimal  sequence  of  basis  functions, and  let 

N  be  an  integer  ^  0  . 

For  any  nontrivial  set  of  vectors  f,g€ff  there  exists  a  constant  C„ 

N 
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(Independent  of  h  }« such  that 


r  i,2N+1/2  ^  t  I  I  h  h  I  I  I 

V  -  -“2Nf''E 

for  h  sufficiently  small. 


Proof. 

The  proof  Is  by  contradiction,  l.e.,  we  assume 

for  some  sequence  h^^  ,  with  h^  0  as  i  ^  . 

If  f  and  g  are  linearly  independent.  Theorem  4.1  of  [6]  then  gives  that 


•/,  \2N+2  _  ,,  ,2N 

span{!|(j}j_Q  ' 


and  this  is  obviously  a  contradiction. 

We  therefore  only  have  to  consider  the  case  when  f  and  g  are  linearly 
dependent,  say  f  «  a-g  ,  g  0  .  As  in  the  proof  of  Theorem  4.1  of  [6]  it  now 
follows  that 


(3) 


d  ,o  d  ,1  £  ,  d  ,  ,2N 

^  Vl  -  “  d7  W«®P^"^d7 


where  i|ij  and  are  as  introduced  in  Lemma  3.1  ot  [6].  Since 


/  d  ■  \2N  , d  .o  d  ,1,N 

8pan{^  '"j-  ^ 


and 


(3)  immediately  leads  to  the  conclus 


'*'1  “  “'*'1 


1  ^ 


Because  of  the  fact  that  i/Zq  and  | 


(A) 


dy  1 


But  according  to  [63  satj| 


k 


i  =■  n 


so  (4)  is  obviously  a  contradlctlc 


For  the  analysis  in  this  paper 


one  concerning  the  true  solution  u 
sionally  reduced  solutions  u|^  • 


m 


i 


T.EMWA  2.2.  Let  u  denote  th«^l 
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h  ^ 

LEMMA  2.3.  Let  'l».(y/h)x.  denote  a  dimensionally  reduced  solution 

j*0  J  J 

of  order  N  .  If  for  some  integer  k  ^  0 


f,g€P((AM  hh 


then 


Xj€0((M=^A)'''^^)  Vj:  0  <  j  <  N 


Proof. 


,N+1 


Let  x€H  denote  the  vector  (xq,..,Xj^)  .  If  is  clear  than  x  is 
the  solution  to 


h<£A^'^^x,A^^^X=’  +  *  <1.Z>  , 

where  C=  ,  D*  “  ^’^l^i-0  ’  given  by 

1  1 

Cij  =  j  ^^(y)^j(y)dy.  d^^  =  |  ^^)dy 


and  r^  -  t{)^(l)g  -  ili^(-l)f  respectively.  (We  have  here  used  also  to 

N 

denote  the  Inner  produce  <x.  ,y,>  in 

1«0  ^ 

Due  to  the  fact  that  A  is  self-adjoint  and  M  bounded  we  get  that 
N+1 


^  *  £  »  i .  e . , 

^  -  h“^C"^r  -  h"^C"^Dto 


X  €[I?(A)] 


and 
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By 

implies 


—1  k  N+1 

successive  application  of  this  equality  it  follows  that  r^€lP((AM  )  )] 

X  €  .  This  finishes  the  proof  of  the  lemma. 

□ 
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3.  A-POSTERIORI  ERROR-ESTIMATION  IN  THE  GENERAL  CASE 

As  already  mentioned  In  the  introduction  one  purpose  of  this  paper  is  to 
derive  a  reliable  technique  for  a-posterlori  estimation  of  the  error  introduced 
by  dimensional  reduction.  The  error  here  is  measured  in  the  energy-norm.  The 
key  Ingredient  of  this  technique  Is  a  socalled  estimator  Est,  which  we  now 
proceed  to  define. 

Let  e  €H^(r-h,h];W)  be  the  solution  to 

a^  A  Me  =  g  _  ^  Mujj  for  y  =  h 

=  f  -  for  y  «  -h  . 


Since  (6)  is  a  Neuman  problem  for  e  it  only  has  a  solution  provided 


h 

<g,x>  -  <f,x>  -  I  bj^<AuJj,x>dy  -  0  Vx^W 
-h 

Because  of  the  equations  defining  u^  it  follows  that  this  is  true  if 

(7)  16span{ii»j}j^Q 

Note.  According  to  Theorem  3.1  and  4.1  of  [6]  condition  (7)  is  in  general 
necessary  and  always  sufficient  to  ensure  that 

I  I  I  Ig  ^  °  ^  ^  °  ' 

Certainly  (7)  is  satisfied  for  any  optimal  sequence  of  basis  functions. 


We  now  define 


(8) 


Est  =  ( 


-h 


The  function  t  is  clearly  not  uniquely  determined,  but  since  only  ^  e  is 
involved,  Est  is  well-defined.  We  shall  now  show  that  the  estimator  Est 
e^ibits  some  very  attractive  properties. 


THEOREM  3.1.  Let  u^  be  the  solution  to  (2)  and  u^  a  dimensionally 
reduced  solution  of  order  N  ^  0  corresponding  to  a  sequence  of  basis  functions 
that  satisfies  (7) .  If  Est  is  as  defined  in  (8)  then 
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1 1  !e  1 

Note.  In  the  terminology  of  [1]  this  theorem  says  that  Est  is  a 
'guaranteed'  upper  estimator. 

Proof 

Clearly 

where  the  sup  is  taken  over  v  €H^( [-h,h] ;H)  nL^( f-h,h] ;P(A^^^) )  .  According  to 
the  definition  of  e  this  is  nothing  but 

h 

supI  I  ^ 

-h  ^ 

and  using  Schwarz's  inequality  we  now  get 

1  i  I  Ig  -  ^  I  l^dy)^''^  =  Est 

-h 

For  use  of  the  estimator  Est  in  actual  computations  it  is  important  that 
it  is  very  close  to  the  real  error  in  a  wide  class  of  situations.  In  the  termin¬ 
ology  of  [1]  this  is  expressed  by  the  requirement  that  Est  be  asymptotically 
exact.  The  following  theorem  contains  a  precise  formulation  of  the  asymptotical 
exactness  for  the  estimator  Est  .  It  is  essential  here  that  the  dimensional 


reduction  is  based  on  an  optimal  sequence  of  basis  functions. 
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THEOREM  3.2.  Let  ^'^j^j=0  optimal  sequence  of  basis  functions.  Let 

u''  be  the  solution  to  (2)  and  u|j  be  the  dimensionally  reduced  solution  of 
order  N  ^  0  .  Assume  furthermore  that  f  and  g  are  elements  of 


P((AM  ") 


fia-X-H 


If  Est  is  as  defined  in  (8)  then 


Est  =  II  I  !g(l+0(h  )) 


Note.  Here  [•]  denotes  the  Integer  part. 


Proof 


Since  f,g€l?(AM  it  follows  from  Lemna  2.3  and  (5)  that  e  in  this  case 
can  be  selected  so  that 


e€H^([-h,hl;P(A)) 


We  now  have 


Est  = 


ft! 

I  „l/2  d  „l/2  d  . 


n 

‘dy|/(  e|  l^dy)- 


n  n 

<  (supl  I  ^  ^  v>dy|/| | iv| i |g) • ( | | |e | | |g/(  j  Sj^j ~  e| |^dy)^^^)  , 


where  the  sup  is  taken  over 


V  ^H^([-h,h] ;«) n  L^([-h,h] ;P(A^^^))  . 


As  In  the  proof  of  Theorem  3.1  this  last  expression  is  equal  to 


h 

Il|»''-u5|||^  •  ||U|||^/(  I 

-h 


The  theorem  will  therefore  be  proven  if  we  show  that  e  can  be  chosen  such  that 


(9) 


11  ^dy  <  Cjjh^ 


1 


It  is  clear  that  by  appropriately  selecting  the  undetermined  constant  of 
t  ,  we  can  obtain 


(10) 


=  Ch 


d 

^^dy 


Aody 


-h 


Now  from  (6) ,  the  definition  of  c  ,  we  have 


^  EyTuT  Ae>dy  » 


h  dy  ’ dy 


^  ^  H'^Aody 


-h 


-h 


h 

*  <g,M"lAE(h)>  -  <f  ,M"^Ae(-h)>  -  |  ^  ^  M'^Aody  - 

-h 

h. 

-  j  bjj<Au^,M"^AE>dy 

-h 


and  Introducing  this  Is,  because  of  Lemma  2.3,  equal  to 

Bj^(M“^A(u*^-ujj)  ,e) 


By  an  application  of  Schwarz's  inequality,  this  expression  can  be  bounded  by 

M  11^  •  nielli^  . 

From  the  fact  that  A  is  strictly  positive-definite  together  with  (10)  it 
follows  that 

lllelllg  <  C(l+h^)(|  aj^ll^  A^^^Ell^dy)^/^ 

-h 

and  since  we  only  need  to  consider  small  h  this  gives 
h 

®hl'^  =  8j^(m"^A(u*'-uJJ),e)  < 

-h 

h 

<  C I  [  iM-^ACu^-uJ)  i  I  I  e(|  %  1  1-^  A^^^e  I  I  ^dy) 

-h 

By  insertion  in  (10)  we  conclude  that 
h 

(11)  I  b^l  i^dy  <  Ch^l  1  |M-1a(u*’-u!J)  I  I  |2 

-h 

For  the  rest  of  this  proof  let  us  assume  that  N  is  even  (the  procedure 
for  N  odd  is  quite  similar,  only  there  are  slight  variations  of  Lemma  2.1  and 
Theorem  3.1  of  [5]  needed  in  this  case). 


It  follows  immediately  from  the  proofs  of  Lemma  2.2  and  2.3  that  M 
and  M~^Au^  are  solutions  to  the  same  problems  as  u''  and  just  with  f 


N 


and  g  replaced  by  AM~^f  and  AH~^g  .  Because  of  Theorem  3.1  in  [6]  and  the 
assumption  that  f,g€0((AM  we  get 


.-1 


I  I  |«-1a(u’’-u5)  1 1  <  C„h 


N+1/2 


and  this  combined  with  (11)  immediately  leads  to 


V  I  1 .1/2  11 2.  ^  U2M+3 


-h 


On  the  other  side  from  Lemma  2.1  it  follows  that 


h  .hi  m2 


From  the  definition  of  Est  and  Theorem  3.1  of  this  paper  we  therefore  get 


0/-^'  i 


-h 


that  is  we  have  finally  proven 


j  bJlA^^^^ll^dy  <  Cj^h^l  a^llM^/^^tll^dy 

-h  -h 


□ 


Since  the  estimator  Est  is  to  be  used  in  actual  computation  it  is  of  utmost 
Importance  that  it  can  be  calculated  very  simply.  The  equations  (6)  and  (8)  that 


19 


Include  the  solution  of  an  0.0. E.  are  therefore  not  well  suited  as  a  formula 
for  the  calculation  of  Est  .  In  the  following  we  shall  show  how  easily  Est 
can  be  calculated  by  means  of  different  formulae.  For  these  formulae  to  be  valid 
it  is  essential  not  only  that  the  dimensional  reduction  Is  based  on  an  optimal 
sequence  of  basis  functions  but  also  that  this  sequence  satisfies  a  special 
orthogonality  condition. 

OD 

Let  optimal  sequence  of  basis  functions  which  is 

orthogonal  in  the  semi  inner-product  f^a  ^  ‘  clearly 

uniquely  determined  modulo  a  constanf  and  a  scalar.) 

oo 

We  define  a  sequence  by 


1  1 

lij^(u)  =■  <|;j^(y)  -  (j  b(t)dt)"^  •  j  b(t)\()j^(t)dt 


(♦jf  j  ^  1  »  3re  therefore  uniquely  determined  modulo  a  scalar.) 


LEMMA  3 . ]  ■  Let  be  a  sequence  as  defined  above.  Let 

h  ? 

”  I  <f'j(y/b)x.  denote  the  corresponding  dimensionally  reduced  solution  of 

j=0  J  J 

00 

order  N  .  If  e  Is  as  defined  In  (6)  then  there  ::xlst  constants 

,  ,00 

such  that 

^  <^l)(y/h)M'^(Co(f+g)  +  c^(f-g))  + 


+  iI>2)(y/h)hM  ^ACj^Xq  for  N  ■  0 
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and 


'•’l«-2^  (y/h)hM~^ACjj^^Xj^  for  N  >  1  . 


Proof 

Since 

imnediately  Chat 


is  an  optimal  sequence  of  basis  functions  it  follows 


tH-2 


e  =  I  (|>  (y/h)E, 
j=0  J  ^ 


for  some  £^€H  0  1  j  £  N+2  ,  l.e. , 


H 


1/2 

For  any  1  £  j  £  N  and  x  CP(A  ) 


1 

<Mej,x>j  a(~  i|ij)^dy  =  ,M^^^x>h 

-1  -h 


a^(^['t'j(y/h)])^dy 


^  [(|.j(y/h)x]>dy 


-h 


where  the  last  equality  is  due  to  the  fact  that  orthogonal  in  the 

semi  inner-product  /^a  ^  ^  from  (6)  we  get  that  the  last 

expression 
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(y/h)x) 


.  h 


and  this  vanishes  because  is  defined  as  a  projection.  We  have  therefore 

proven  that 


1  /2 

<Mej,x>  =  0  Vl  1  j  1  N,  )  , 


or 


(12) 


Cj  =  0  Vl  <  j  <  N 


For  j  =  N+1  we  get  as  before 


<Mej^l.x> 


-1 


-h 


and  this  is  according  to  (6) 


=  h<g,<frjj^^(l)x>  -  h<f,4ij^j(-l)x>  -  hj  bj^<Au5j,(t.j^^j(y/h)x>dy 


-h 


Jn  this  identity  we  also  used  the  orthogonality  of  the  "fj'®  •  Concerning  the 
last  term 

h  N  ^ 

^  j  I  <AXj,x>  I  b(y)(j>j(y)4ij^j^(y)dy  . 

-h  -1 

This  vanishes  for  N  -  0  ,  because  /^b(y)*. (y)dy  «  0  .  and  therefore  establishes 

-1 

the  formula  for  in  the  particular  case  N  -  0  . 

In  the  following  we  are  left  to  consider  N  ^  1  .  Let  j  be  any  Integer 


0  £  j  <  N-1,  then 
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values  in  the  case  of  constant  coefficients  a  and  b  . 


Example  3.1 

If  the  functions  a  and  b  are  both  constant  it  immediately  follows  that 

00 

every  optimal  sequence  of  basis  functions  satisfies 


span{ii/j  =  all  polynomials  of  degree  N 


for  any  N  ^  0  ,  and  vice  versa. 

The  specific  optimal  sequence  of  basis  functions 
3.1  and  Theorem  3.3  is  now  (modulo  a  scalar) 


used  in  Lemma 


*  y 

<tj(y)  =  /^^j_j^(t)dt  ,  j  1  2  , 

where  denotes  the  Legendre  polynomial  of  degree  k  . 

Because  of  the  alternating  even  and  odd  polynomials,  it  immeidately  follows 

that 


Cr  =  =  0  Vk  1  0 

Simple  algebraic  manipulations  with  the  Legendre  polynomials  now  give 

Cq  ■  l/2a 
Cj^  *  b/a  and 

c^^  =  b/((2k-l)(2k-3)a)  for  k  >  2 
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From  the  definition  of  Cj^  immediately  follows  that 


1  H 

Ck  *  W  -^k  •  i*"'* 

Cq  =  [1/2  a  ] 

=  [2b^/3a]^'^^  and 

9  ?  2 

Cj^  =  [2b  /(2k+l)(2k-l)‘^  2k-3)  a)]  for  k^2 


We  shall  now  show  how  one  can  derive  another  set  of  formulae  for  the 
a-posteriori  error  estimator.  As  it  will  turn  out  these  formulae  are  much  better 
suited  for  practical  applications.  For  the  case  N  =  0  it  is  clear  that 
Axq  =  Ch~^(f-g)  .  In  the  following  we  therefore  only  consider  N  ^  1  . 


LEMMA  3.2.  Let  ^S^k=0  ^  ^  previous 

11 

theorem.  Then  there  exist  constants  1  <  i  ,  j  2  such  that 

N  —  — 


+  h-loJ^MCf-a^  ^  Mujjc-h))  and 


+  h'^D^HCf-Sj^  i  MnJjC-h))  for  K  >  1 


Proof. 


From  (5)  and  Lemma  3.1  it  follows  that 
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■  ^dy  ®  dy  *<’n+2^  ^‘^N+l’S^ 


If  we  integrace  the  right  hand  side  of  (6)  by  parts,  set  v  -  4iQ(y/h)*x  and 
apply  the  identity  (16),  the  result  is 


hdjjSl  A(Cj^Xjj_j^+CjjXjj)  +  hdj^^  g  f  -  a^  ^“wl-h 


X 

‘^N  =  J  d7  ^  ^  ^N+i‘*y 


jl2  d  d  .  , 

‘^N  "1  dy  ®  dy  ’’’n+z'^^ 


Performing  the  similar  procedure  with  v  =  (<)j^(y/h)'x  instead  we  get 


hdjj  M  ^A(CjjXjj_j^+Cj^Xjj)  +  hdjj  M  ACj^^^x^ 


J. 

^21  f  d  d  .  .  , 

‘^N  *  J  ^  ^  ^  Vi'^i'^y 

-1 

^22  f  d  d  ^  .  . 

‘^N  ’  J  ^  ®  d7  ♦N+2<^i‘^y 
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The  lenuna  therefore  immediately  follows  if  the  matrix 


In  “  ^V^i,j=l 


is  always  invertible.  The  proof  of  the  invertibility  follows  by  contradiction. 
Assume  ^  is  not  invertible  for  some  N  .  This  implies  that  some  nontrivial 
linear  combination  of  the  columns  of  d„  vanishes.  In  terms  of  the  functions 
and  says  that  some  nontrivial  linear  combination 

®2‘^IW-2  such  that 

1 

(17)  I  ^  ^  ^  <^1<^I«-1'^®2W^  V  dy  =  0 

-1 

for  V  *  <j)^  and  v  =  ({i^^  .  By  performing  an  integration  by  parts  we  easily  see 
that  this  Identity  must  also  hold  for  v=ij(^  2£j£N.  Since  the  sequence 
is  an  optimal  sequence  of  basis  functions  we  know  that 


dy  ®  dy  ^®l‘'’w-l'^®2*N+2^ 


tl 

for  some  set  of  constants  (®j^j=0  '  Combining  this  with  (17)  we  conclude  that 


^  ^  <®lVl-"®2 V2>  * 


Due 


to  the  fact  that  ®2^N+2^^^  *  ®  ^  i  1  this  implies 


®l'^N+l*‘^  ®2^N+2  “  ® 


'wmt  > 


mmm 
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which  obviously  contradicts  the  fact  that  this  is  a  nontrivial  linear  combination. 


The  difference  between  the  formulae  given  in  Theorem  3.3  and  those  that 
are  based  on  Lemma  3.2  is  that  while  the  first  include  elements  of  the  form 
AXj  the  latter  are  expressed  solely  in  terms  of  f,g  and  the  Xj's  •  In 
practical  applications  we  seldomly  know  the  exact  values  of  the  *j'®  •  I*'” 

stead  we  compute  some  approximate  values  x^  ,  e*g>  using  a  finite  element 
method.  The  error  introduced  by  using  approximate  values  derived,  say  from 
finite  elements,  in  the  expressions  of  Lemma  3.2  can  be  neglected.  The  reason 
is  that  the  difference  between  Xj  and  in  the  H  norm  (viz  L  )  is 

normally  very  small.  The  problem  with  the  expressions  of  Theorem  3.3  in  this 
context  is  that  in  general  AXj  is  not  at  all  defined. 

ii 

As  before  we  now  give  the  values  of  the  constants  in  the  case  where 

a  and  b  are  both  constants. 


Example  3.2 

Assume  a  and  b  are  constants.  In  this  case  we  already  know  that  0 

Lemma  3.2  therefore  reduces  to  formulae  for  ACj^Xjj_j^  and  AC^ ^  ^^x^  ,  N  ^  1  . 

It  is  easily  seen  that  these  are 


ACj,x^_^  =  h‘^2a(2W-l))'^^2lM(g-aj^  ^  MuJ(h))  + 

+  (-l)Vf-a^lpMuJ(-h))] 


ACn+1*n  “  h■^(2a(2^^f3))  Mu^(h))  + 
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4.  IMPROVED  ERROR  ESTIMATION  —  A  SPECIFIC  EXAMPLE 

As  mentioned  before  It  Is  very  Important  that  we  are  able  to  estimate  the 
error  accurately.  Theorem  3.2  shows  that  our  estimator  Est  does  exactly  that 
provided  the  data  is  sufficiently  regular  and  h  is  not  too  large.  In  this 
section  we  shall  address  the  problem  of  how  to  detect  if  the  estimator  Est 
is  too  conservative,  due  to  singularities  in  the  data  or  large  h  ,  and  what 
can  be  done  to  correct  It.  For  simplicity  we  consider  the  model  problem 

=  0  in  »  ]0,1[  X  ]-h,h[ 

=  g  for  y  =  h 

3u^  f  I, 

•  -g  y  “  -h 

u^  -  0  for  X  “  0,1  , 


where  g  is  an  element  of  L  ([0,1])  . 

00 

Let  sequence  of  polynomials  introduced  in  Example  3.1. 

The  dimensionally  reduced  solution  of  order  N  ,  “jj  *  f’a®  Ihe  form 


.  N 

UvCx.y)  -  I  <[>.(y/h)v  (x) 
j>0  J  ^ 


o  ^ 

where  Vj€H  ([0,1]),  0  ^  j  ^  N  .  Let  e^^  denote  the  exact  error,  i.e., 
Bjj  ■  u^  -  u^  .  Since  everything  is  even  in  y  all  the  terms  in  the  dimen¬ 
sionally  reduced  solutions  corresponding  to  odd  indices  vanish.  From  here 
on  we  only  consider  dimensionally  reduced  solutions  of  even  order  2N  . 

Let  us  start  by  giving  a  table  that  shows  the  relative  error 


(»1 1  |e2jjl  I  Ig/ 1  1  U^l  1  Ig)  the  efficiency  index  of  the  estimator  Est  (Eff  = 

1  I  |e2jj I  I  Ig/Est)  in  the  case  g(x)  =  it/4  ,  for  N  =  0  and  N  ■  1  respectively, 
and  for  different  values  of  h  . 


Table  4.1 

N  =  0:  N  -  1: 


h 

Rel. Error 

Eff. 

h 

Rel. Error 

Eff. 

1/2 

0.67 

0.88 

1/2 

0.121 

0.69 

1/4 

0.43 

0.94 

1/4 

0.051 

0.68 

1/8 

0.24 

0.97 

1/8 

0.019 

0.68 

1/16 

0.12 

0.99 

1/16 

0.007 

0.68 

1/32 

0.06 

0.99 

1/32 

0.002 

0.68 

From  the  table  it  is  evident  that  the  efficiency-index  approaches  1,  as 
h  0  ,  for  N  =  0  ,  but  that  this  is  not  so  for  N  =  1  .  The  numbers  therefore 
clearly  show  that  some  smoothness  condition,  as  in  Theorem  3.2,  is  essential  in 
order  to  ensure  that  this  index  converges  to  1  for  h  ->•  0  .  For  most  practical 
applications  though,  an  efficiency-index  of  0.7  is  completely  satisfactory  and 
no  corrections  to  Est  are  needed.  It  is  also  important  to  note  that  h  “  1/2 
corresponds  to  a  square  and  that  Est  still  gives  a  very  reliable  estimate  for 
the  error.  The  next  table  lists  the  efficiency-index,  also  in  the  case  where 
g(x)  «  Tr/4  ,  but  for  N  =  2  and  3  and  3  different  values  of  h  . 

By  a  comparison  of  Tables  4.1  and  4.2,  it  is  seen  that  the  efficiency  sig¬ 
nificantly  decreases  as  we  Include  more  and  more  polynomials. 
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Table  4.2 


N  =  2: 

N 

=  3: 

h 

Rel. Error 

Eff . 

h 

Rel. Error 

Eff. 

1/2 

0.045 

0.53 

1/2 

0.022 

0.43 

1/4 

0.019 

0.52 

1/4 

0.008 

0.38 

1/8 

0.008 

0.55 

1/8 

0.003 

0.41 

In  the  following  we  shall 

take  3 

i  closer  look  at 

the 

derivation 

of  the 

estimator  Est 

for  the  purpose 

of  suggesting  corrections 

that  can  increase  the 

efficiency  to 

any  desired  level 

.  We 

shall  only  work 

out 

the  details 

of  a  first 

correction.  ^ 

The  exact  error  e2jj  Is  the  solution  to  the  boundary  value  problem 

^3y^  ^®2N  ’^2N  S 

^®2N 

~ir  °  '’2N  y  *  ^ 

^®2N 

- P2J,  for  y  =  -h 

e2N  =  °  • 

with 

ri'3'k2..3.2,  h,  - 

^2N  =  <37^  ]u2N(’^.y> 


P2n  “ 


J-  .h 

3y  2N 


(x,h) 


and 


In  terms  of  complimentary  energy  the  norm  of  e2jj  ,  |  !  |62nI  ^ 
characterized  by 

1  h 

ll|e,  lllp=  min  [  [  (8^{x,y)+t^(x,y)N  dydx  , 

where 


(8,t)€M  Iff 

(s.t)  , 

T~  3  +  -5—  t  *  r~.,  , 

3x  3y  2N 

~  ^2N  y  =  h  and 

t  =  -P2jq  for  y  =  -h 

If  we  define  t^  by 

3y  *^o  “  *’2N 

t  *  P^vi  for  y  *  h  and 

O  ZN 

^  *  -P2N  y  * 

Eat  ■ 

functional. 


It  Is  clear  that  (0,t^)€M  .  On  the  other  hand  It  Is  also  clear  that 
/*'  t^(x,y)dydx)^^^  ,  l.e.,  Est  Is  simply  a  particular  value  of  a 

A  «-  O 


3. 


the  minimal  value  of  which  is  the  exact  norm  of  the  error.  One  way  to  improve 
the  estimatior  Est  is  therefore  to  take  the  minimum  over  more  than  just  the 
single  function  t^  .  This  should  not  be  exagerated  since  we  also  have  to  keep 
the  formulae  simple.  Define 


s(x,y) 


X 

(r2j,(z,y)  -  r2jj(2)(^  C)(y/h))dz 

O 


and 


t(x,y)  =  hr^jjCx)^  (y/h)  , 


where  r2jj(x)  denotes  the  function 

h 

-h 

and  c  is  an  arbitrary  element  of  •  It  is  obvious  that 


If  furthermore  ^  satisfies  ^(1)  =  1  ,  C(-l)  =  -1  it  also  follows  that 


t  =  P2jg  for  y  =  h  ,  and  t  =  -P2j,  for  y  =  -h 


From  Lemma  3.1  we  get  the  identity 

t2j}(x,y)  “  (-^  » 

where  denotes  the  Legendre  polynomial  of  degree  2N+1  ,  so  it  Immediately 


follows  chat 


1  h 


. 


(s^(x,y)  +  t^(x,y))dydx 


0  -h 


1  X 


(h- 


^  (y)dy 


(r2jj(x))  dx  +  h 


-1 


( 

;  j 
0  0 


r2jj)^dx) 


1/2 


It  should  be  noticed  that  the  estimator  Est  Is  obtained  from  this  with  the 
choice  ?(y)  =  *-210+1  •  If  we  define 


A(h) 


B(h) 


1 

h^j(r2jj(x))^dx  and 
0 

1  X 

h  |(|r2j,)  ^dx  ,  then 
0  0 


the  previous  expression  can  be  written  as 


(18) 


1 

2 

(  C  (y)dyA(h)  + 
-1 


r  a  9  1/2 

J  V'^2N+i(y>-^<y»> 

-1 


It  Is  easy  to  see  chat  the  minimum  of  this  expression  over  Q  approaches 
1 

(  “Est  as  h  -*■  0 

-1 

provided  A(h)/B(h)  -►  0  for  h  -+  0  .  In  the  case  that  A(h)/B(h)  -♦  *  for 
h  -»■  0  Che  minimum  of  the  expression  (18)  approaches 


(4A(r.)B(h))"^^  as  h  -►  0  . 


Based  on  these  asymptotics  we  Introducs^ 


.■{ 


Est  if  A(h), 


(4A(h)B(h)) 


where  t  is  some  specified  constant, 
tell  us  whether  we  shall  use  the  value 
Since  we  know  that  hr2jj(x)  = 
expressed  in  terms  of  P2jj  .  It  folloulf 
exactly  if  Is  smooth.  Since  the' 


the  interval  [0,1]  we  get  that  the 
tion  of  the  error  in  parts  of  [0,1] 


the  interval  [0,1]  into  K  disjoint' 


Est,  =  (  I 


where  Estj^(Ij)  refers  to  the  estimatasi* 
The  following  table  shows  the  efflciettCJ^ 
I  I  ie2j^|  I  |j./Est2)  in  the  case  g(x)  -  »/| 
values  of  h  .  Tne  interval  [0,1] 
length,  l.e.  K  =  8  ,  and  the  constant  f 


chosen  to  be  10 


Although  the  efficiency-index  is 
improvement  over  Table  4.2.  Again  nOiat  j 
of  h  .  If  additional  accuracy  of  tlM 
obviously  be  obtained  by  an  extension  ' 


Whether  such  additional  corrections 
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5.  S(»1E  REMARKS  ON  AN  ADAPTIVE  STRATEGY 

As  already  mentioned  in  the  Introduction,  the  goal  of  this  paper  Is  not 
only  to  derive  reliable  estimates  for  the  error,  but  also  to  use  these  estimates 
as  tools  in  an  automatic  selection  of  the  right  order  dimensionally  reduced 
solution  for  a  given  problem. 


First  let  us  introduce  a  slight  generalization  of  the  concept  of  dimen¬ 
sionally  reduced  solution.  Instead  of  projecting  onto  the  space 
N  ,  N 

{  I  <('^(y/h)x  Jx,  €2?(a  '  )  ,  j  =  0,..,N}  we  project  onto  {  J  i/;.  (y/h)x.  |x.  €  K  , 

J  J  J  j=0  ^  ^  ^ 

j=0,..,N}  ,  where  ^^j^j-0  ^  family  of  closed  subspaces  of  V(A^^^) 

To  see  the  importance  of  this  generalization  and  describe  the  ideas  behind 

the  self-adaptive  strategy  we  shall  consider  the  case  that  A  is  a  differential 

operator  on  some  domain  .  Let  be  divided  into  k  disjoint  subdomains 

,  1  £  i  i  k  »  and  let  l£i^k,  be  k  nonnegative  integers.  Set 

1/2 

K.  =  {u€0(A  )  (u(x)  =  0  for  x€  U  R  }  ,  the  extended  concept  of  dimensionally 

^  N^<j 

-  N 

reduced  solutions  with  this  family  {K  } ,  „  ,  N  =  max{N.}  ,  is  one  that  permits 

1  j=0  ’  i  1  ’ 

different  order  of  the  dimensionally  reduced  solution  in  different  parts  of  the 
domain  R  .  This  is  extremely  important  for  practical  applications,  where  a 
low  order  dimensionally  reduced  solution  may  very  well  be  satisfactory  in  the 
Interior  of  the  domain  and  away  from  singularities  in  the  loads  and  at  the  same 
time  a  high  order  solution  is  required  near  the  boundary  or  near  singularities. 

As  a  total  estimator  for  the  error  let  us  use  an  expression  of  the  form 


k  2  ^''2 

(  I  [n,(N.)r) 

i=l 


(19) 
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where  refers  to  some  estimator  on  the  domain  with  respect  to 

dimensional  reduction  of  order  could  for  example  be  the  estimator 

Est  or  the  corrected  estimator  Est2  of  the  previous  section.) 

If 

To  set  a  goal  for  the  'best'  distribution  of  the  orders 
a  dimensionally  reduced  solution  we  need  the  concept  of  cost.  Let  us  assume 
that  the  cost  of  (solving)  the  dimensionally  reduced  problem  with  orders 
given  by 

k 

[  (6N  +l)“m(fi  )  , 

1=1 

where  a  and  0  are  two  positive  constants  and  m(n^)  is  some  measure  of  . 

As  a  'best'  distribution  of  the  orders  for  a  dimensionally  reduced  solution 
we  define  one  which  for  a  given  cost  minimizes  the  energy  norm  of  the  error. 

(We  could  also  have  defined  a  'best'  distribution  as  one  that  for  a  given  value 
of  the  energy  norm  of  the  error  minimizes  the  cost.  Which  of  these  two  defini¬ 
tions  we  take  makes  no  difference  in  the  strategy  we  propose.) 

The  following  is  very  heuristic  in  nature  and  by  no  means  an  exact  verifi¬ 
cation  that  the  strategy  works.  Let  us  use  the  expression  (19)  as  if  it  were 
the  exact  norm  of  the  error.  Secondly  let  us  assume  this  expression  to  be  defined 
for  all  positive  values  of  the  N^'s  ,  and  not  only  Integers.  By  introduction  of 
Lagragean  multipliers  it  is  easily  seen  that  a  'best'  distribution  of  the  N^'s 
has  to  satisfy 


3  C  (Independent  of  i  )  such  that 

(BN^+l)“"^m(nj) 


(20) 


C  ,  Vi 


AO 


In  practice  we  only  have  the  values  of  the  Integer  points  and  a 

discrete  equivalent  of  (20)  Is  then 


(21) 


[n^(N^+l)]^  -  [n^(Np]^ 
(6N^+l)“"^m(Q^) 


C  , 


Vi 


We  shall  also  assume  that  n(N^+l)  is  significantly  smaller  than  n(N^)  so 
that  Instead  of  (21)  we  get 


(22) 


[n^(N^)]^ 

(eN^+l)“"^m(Q^) 


Vi  . 


The  strategy  we  propose  is  one  that  alms  at  equilibrating  the  left  hand 

sides  of  (22).  We  do  this  in  a  way  similar  to  the  adaptive  finite  element  solver 

F.E.A.R.S.  (cf.  [1])  .  Let  us  assume  that  we  have  arrived  at  a  distribution 
o  k 

and  that  the  estimate  for  the  error  is  unacceptably  large.  Our  strategy 
is  simply  to  find  j  such  that 

[n.(<)]2 

6  =  - i — i- - 

^  (6N°+1)®"  m(fij) 

is  maximal,  and  then  Increase  Nj  by  1  .  In  the  next  section  we  shall  see 
how  well  this  performs  in  a  practical  example. 
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6.  A  NUMERICAL  EXAMPLE 

Consider  the  same  problem  as  in  section  4,  namely 


3.2  .  ,  3.2I  h 


0  in  “  ]0,l[x]-h,h[ 


for  y  =  h 


-g  for  y  =  -h 


h 

u 


0  for  X 


0,1 


Let  [0,1]  be  divided  into  the  four  subintervals  =  [(i-l)/4,  i/4] 

1  1  ^  4  .  A  dimensionally  reduced  solution  can  nov  have  different  order  in 

the  different  intervals  ,  1  £  1  4  .  As  basis  functions  we  choose  the  poly¬ 

nomials  Introduced  in  Example  3.1. 

The  equations  that  define  the  dimensionally  reduced  solutions  are  solved  by 
introducing  a  Finite  Element  discretization  In  the  x-dlrectlon.  Piecewise  linear 
functions  on  a  regular  mesh  are  used  as  test  and  trial  functions  for  the  Finite 

Element  Method.  Since  we  want  to  illustrate  the  behavior  of  the  dimensional  re¬ 

duction,  and  are  here  not  interested  in  any  contribution  from  the  x-dlscretization 

_9 

we  choose  a  very  fine  grid  of  meshslze  =  2  .  The  involved  linear  equations  are 

solved  by  a  Cholesky  decomposition  combined  with  iterative  refinement.  In  the 
computations  that  we  present  here  g(x)  is  chosen  >  7r/4  .  Since  this  choice  of 
boundary  data  makes  the  problem  symmetric  In  the  line  x  ^  we  only  need  con¬ 
sider  X  in  the  Interval  L0,H1  •  Let  2Nj^  ,  1  £  i  ^  2  ,  denote  the  order  of  the 

dimensionally  reduced  solution  In  1^  .  The  following  table  shows  the  error  on 

the  whole  Interval  [0,1]  (■  ||le  _t||_)  and  the  work  (here  defined  by  >i(N,+l)  + 

211  ^ 

3  _ 

^(M^+l)  )  as  a  function  of  the  pair  N  >  (N^,N2)  for  two  different  values  of  h 
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Table 


6.1 


2N 


1 

0 


.2N, 


2N 


1 

0 


0 

2 

4 

6 

0.3974 

0.3257 

0.3168 

0.3160 

1 

4.5 

14 

32.5 

0.3662 

0.0730 

0.0721 

0.0720 

4.5 

8 

17.5 

36 

0.3517 

0.0383 

0.0270 

0.0270 

14 

17.5 

27 

45.5 

0.3499 

0.0292 

0.0134 

0.0134 

32.5 

36 

45.5 

64 

0 

2 

4 

6 

0.3014 

0.2296 

0.2270 

0.2267 

1 

4.5 

14 

32.5 

0.2540 

0.0359 

0.0359 

0.0359 

4.5 

8 

17.5 

36 

0.2484 

0.0134 

0.0130 

0.0130 

14 

17.5 

27 

45.5 

0.2479 

0.0063 

0.0055 

0.0055 

32.5 

36 

45.5 

64 

error 

work 
h  = 


error 

work 

h  »  V 


Based  on  the  numbers  In  this  table  we  can  now  find  the  entries  with  the  pro¬ 


perty  that  the  error  is  smaller  than  any  error  obtained  with  the  same  or  less 
work.  These  entries  are  marked  in  the  following  table. 


The  tables  6.1  and  6.2  clearly  illustrate  the  advantage  of  a  non-unlfonn 
distribution  of  the  pol)monlals.  It  is  easy  to  see  that  the  true  solution  u 
in  the  limit  as  h  -*•  0  hag  a  parabolic  behavior  in  the  y  direction,  also 
for  X  in  the  middle  of  the  interval  C0,l].  This  is  reflected  in  the  fact 
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that  the  pair  (0,2)  is  slightly  better  than  (2,0),  it  also  explains  the  signifi¬ 
cant  decrease  in  the  error  obtained  by  choosing  the  pair  (2,2).  For  the  higher 
order  polynomials  there  is  a  clear  tendency  towards  concentration  near  the  boun¬ 
dary  X  *  0  (and  x  =  1)  in  the  entries  marked  in  Table  6.2  .  This  concentration 
is  more  visible  the  smaller  h  is  ;  for  h  »  *5  the  pair  (6,2)  is  not  as  good  as 
(4,4)  but  for  h  =  1/4  the  error  obtained  by  (6,2)  is  less  than  half  the  error 

by  (4,4)  with  only  a  slight  Increase  in  the  work. 

We  now  want  to  test  the  adaptive  strategy  outlined  in  section  5  on  this  exam¬ 
ple.  We  consider  the  case  h  =  ^  ,  where  non-uniformity  in  the  distribution  of  the 

polynomials  is  most  advantageous.  As  an  estimator  we  use  Est2  of  section  4, 

2 

with  the  constant  t  set  to  10  and  each  interval  I  ;  divided  into  2  subin¬ 
tervals  of  length  1/8.  The  following  table  shows  the  efficiency  of  Est2  (Eff2 

1  1  |e2ul  1  I  E/Est2)  as  a  function  of  the  pair  2N  =  (2Nj^,  2N2)  . 

V2N 

Table  6.3  2nX  0  2  4  6 


0.94 

— 

0.99 

1.00 

1.00 

1.06 

0.68 

0.68 

0.68 

1.09 

0.60 

0.59 

0.59 

1.09 

0.57 

_ 

0.53 

0.53 

h  »  it 


It  is  evident  from  Table  6.3  that  Est2  provides  a  reliable  estimate  for  the 
error  even  in  the  case  of  variable  order.  We  also  note  that  Est2  Is  not  nec¬ 
essarily  an  upper  bound  for  the  error,  when  the  orders  of  the  polynomials  are 
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allowed  to  vary.  Steps  could  be  taken  to  correct  this,  but  on  the  other  hand  com¬ 
putational  experience  shows  that  this  effect  is  insignificant,  and  that  Est2  is 
very  close  to  an  upper  bound  in  most  cases. 

Let  us  start  with  an  initial  distribution  for  the  orders  of  the  polynomials 
given  by 


(2Nj^,2N2)  ■=  (0,0) 

Based  on  the  present  formula  for  the  work  and  the  error  estimate  we  now  compute 
,  j  =  1,2,  as  in  section  5.  The  result  is 

*  ^2  “ 

We  can  therefore  proceed  to  both  (0,2)  and  (2,0).  According  to  Table  6.1 
(0,2)  is  only  slightly  better  than  (2,0),  so  this  apparent  "failure"  of  our  strat¬ 
egy  is  of  very  little  significance. 

For  the  pair  (0,2)  we  compute 

6^  -  0.10  ,  ^2  “ 


and  for  the  pair  (2,0) 


6^  -  0.31  X  lO'^  ,  62  “ 

In  both  of  these  two  cases  we  are  told  to  proceed  to  the  distribution  given  by 


(2,2) 
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For  this  pair  we  get 


0.14  X  10 


-2 


and  6, 


0.59  X  10 


-6 


l.e.,  if  ve  want  higher  accuracy  with  dimensional  reduction  our  strategy  selects 
the  pair 


(4,2) 


In  this  case 


=  0.14  X  10  ^  and  62  =  0.16  x  10  ^  , 

so  that  additional  requirements  to  the  accuracy  will  lead  us  to  the  dlstributio  . 

(6.2) 

The  path  that  our  strategy  goes  through  can  schematically  be  represent!,  j  as 

.  (0.2)  , 

(0,0)  ^  (2,2)  (4,2)  (6,2)  , 

(2,0)  ^ 

and  based  on  the  tables  6.1  and  6.2  this  Is  clearly  seen  to  be  a  very  good  choice. 
The  strategy  has  been  tried  In  a  variety  of  other  situations  and  has  consistently 
been  very  effective.  It  has  also  been  tried  with  different  measures  for  t'l.e  work. 
Here  It  should  be  noted  that  by  changing  the  measure  of  the  work  we  may  entirely 
change  the  "best"  distributions  for  the  polynomials,  but  the  strategy  detects  that 
easily . 
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7 .  CONCLUSIONS 

In  the  following  we  list  some  conclusions  concerning  the  approach 
of  dimensional  reduction  developed  In  a  series  of  three  papers  ([63,  [7],  and  the 
present) . 

a)  It  Is  common  In  engineering  to  distinguish  between  structures  with  large 
and  small  thickness  (see  e.g.  [3]).  The  approach  presented  here  entirely 
avoids  this  somewhat  artificial  categorization. 

b)  This  approach  gives  in  an  optimal  and  adaptive  way  the  advantages  of  asymp¬ 
totic  expansion  (when  the  thickness  is  small)  and  the  effectivity  of  the 
spectral  or  p-verslon  methods  (when  the  thickness  is  not  small,  or  strong 
singularities  are  present).  It  has  been  shown  that  these  two  requirements 
uniquely  characterize  the  approach. 

c)  Reliable  aposterlorl  error  estimates  can  be  obtained  for  this  approach,  and 
they  lead  immediately  to  an  effective  adaptive  strategy. 

d)  The  approach  Is  numerically  very  robust  and  works  well  independent  of  the 
thickness  and  regularity  of  input  data. 

e)  The  underlying  mathematical  theory  and  numerical  experiments  show  the  direc¬ 
tion  for  various  generalizations.  These  shall  be  dealt  with  elsewhere. 


i 
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